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Abstract Globally optimal solution describing a phase conjugated field of Raman scatter- 
ing on the resonant B «— X transition of iodine I2 is studied. Maximum optical 
coherence is found as a top eigenvalue problem. A reversibility theorem has 
been stated. This provides sufficient conditions for a tightly localized waveform 
and molecular hologram to exist. A noisy picosecond pulse has been computed 
to show how femtosecond polarization is regained at target time. 

1. Introduction 

Molecular wave packet engineering has attracted much attention in the works 
collected under the rubric of "quantum control" 1_14 . Current world-wide ef- 
forts in the problematic have been mounted to develop efficient methods for 
breaking selected molecular bonds 11-14 or to harness specific molecular states 
for optical processing devices and spectroscopic uses. The rapid progress did 
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not take a long time, because it was prepared by the enormous lore in photo- 
molecular spectroscopy accumulated since Lord Rayleigh epoch and the be- 
ginning of quantum mechanics. From the other hand, optimal control theory 
comprising dynamic programming and modern variational calculus has been 
the subject of mathematical studies enabling to propose a theoretical apparatus 
to the quantum control. 

Of prime importance were the minimum quantum uncertainty states intro- 
duced by Schrodinger. For decades they were meant of only as "Gedankenex- 
perimente". Up-to-date femtosecond technique has made it possible to observe 
both classically confined states of Rydberg atoms 4 and space-localized wave 
packets in molecules 6 ' 9 . The title of this paper obliges us to restrict ourselves 
by the latter. The quantum control theory of ultrafast events close to dissoci- 
ation limit of diatomic molecules is our major concern. The challenge here is 
to find optimal laser excitation of picosecond scale causing femtosecond radi- 
ation of an optically thin sample. 

We shall deal with the subject regarding the iodine I2 molecule as our test 
example. Diatomic iodine in gas or condensed state has become the reference 
standard 15 and ideal benchmark 16,17 in modelling the wave packet evolution. 
Our task is greatly facilitated by the considerable volume of researches, in 
which the quantun control of molecular motion has been exploited theoretically 
5 ' 7,9 . The following experiments 6 ' 8,9 with the laser-induced fluorescence (LIF) 
have supported the idea of wave packet localization inside a molecule. These 
works have again emphasized the link of the phase modulated photoexcitation 
and molecular vibrations, which was broadly known in Raman spectroscopy 
from 1920s 18 . 

The goal yet achieved in the dynamic quantum control 9 includes localiza- 
tion of the vibrational wave packet at the attractive side of molecular potential. 
This scheme was referred to as the molecular reflectron 5 . Bulk literature has 
been devoted to focusing the matter states 3-9 . The right posed theoretical 
limits 10,11 indicate an existence of femtosecond coherence. 

Our approach to the problem of a drastical shortening optical transient of a 
single molecule is started from a crucial relation between a time-reverse molec- 
ular dynamics and phase conjugate resonance scattering field. A considerable 
interest presents a justification of the relation, which, despite its generality, we 
have not been able to find in the literature. There are two aspects: designing a 
right chosen objective and tailormaking optimal optical pulses. To understand 
how they appear, it is worth noting, that the rapid improving of femtosecond 
techniques is based on compression of chirped light pulses provided by the 
well-established wave guides in optics. One can borrow the key element of the 
pulse compression physics to apply it to squeezing an optical coherence inside 
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molecular space. Namely, the Franck-Condon region must play a role both dis- 
persive wave guide and frequency modulator owing to molecular vibrations. 

The vibrational wave packet is expected to be focused at will on inner 
Franck-Condon region, since the Franck-Condon factor is at its maximum 
at the steep repulsive curve. Simultaneously the momentum variance Ap « 
2Mvq will be maximal at this point with vq being the wave packet velocity. 
The uncertainty principle ApAq ~ % guarantees a tight localization Aq of 
the rebound wave packet. Thus, its overlap with the ground state, having the 
variance Aqo 3> Aq, lasts just for recoil lifetime (~ Aqo/vo) and gives rise to 
ultrashort coherent transient. This picture appeals to the "billard ball" model, 
19 which still remains to be extended to involve spreading wave packets on 
molecular curves. As example we shall look at the reflectron scheme. 

The specific iodine reflectron 5 ' 6 ' 9 works at certain frequency above the 
ground state X. A tailored electric field excites the B state until the wave packet 
begins to concentrate near an outer turning point far removed from the loca- 
tion of the original Franck-Condon transitions. The closer the excitations to the 
dissociation limit, the longer a delay before the wave packet will be reflected 
from the outer curve and moves back to inner repulsive core. From here the 
wave packet recoils and may radiate photons of resonant frequency to the B-X 
transition, so that, an optical coherence begins to appear with a some delay. 
This "dark" period might cause a fluorescence to spark at a chosen time. 

Having the name drawn from the electronic prototype, the molecular re- 
flectron is the scheme to create vibrational wave packet and to focus it on a 
desired material target. Our prime interest is a radiation process, in which 
the wave packet is periodically returning into resonance Franck-Condon re- 
gion and recoiling from it. Hence, a molecular "magnetron" is the best-fitting 
term to our design, which underlines a parallel between optical-molecular and 
radio-electronic phenomena. In fact, the vibrational wave packet spreading as- 
sisted by molecular anharmonicity could be negated by the tailored laser pulse, 
because an interplay between its modulation and quantum dispersion squeezes 
the wave packet, as does a properly chirped radio pulse propagating along a 
dispersive delay line. 

Molecular states promoted nearby their dissociation limit can be localized 
more tightly than the ground state by itself. For the excited levels are quasi- 
classical in their nature. Thus, we can treat the recoil lifetime r rec as a kine- 
matical overlap of the state X and B presented schematically in Figure 1. Here, 
the resonance Franck-Condon transitions are drawn by the vertical arrows; the 
horizontal arrows designate the wave packet motion. The wave packet velocity 
v and recoil time r rec are found from the potential functions (molecular curves 
shown in Figure 1) as 
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Ro+AR 

v(r) = [2(e-U B (r))/M] 1/ \ r rec ps | dr/t;(r). 

-Ro 

The distance i? = is the equilibrium molecular position in the ground 
X state of I2. The AR is a maximum variance for the overlapping states, 
where e is the energy of excitation (equal to an optical frequency) and M 
being the reduced mass of I2. The variance AR^H/MQ) 1 / 2 ps 0.05 ,4 (with 
w 214cm _1 ) for the ground state is more than it for focused B wave packet 
on the inner molecular wall. Furthermore, if its variance is disregarded, the 
recoil lifetime is limited below by the magnitude r rec ps 15 fs. The state 
B can be populated at the energy e ps 19050cm -1 (~ 525 ran) for a more 
long time than the recoil time r rec . Setting the gate of excitation by the typical 
vibrational period T ps 546 fs at e, one can measure the squeezing by the 
ratio T/r rec which is expected about 30 or more. As it will be clear later, 
even subtle details of controlling the atomic motion inside diatomic molecules 
can be understood by means of the classical characteristics and semiclassical 
distribution functions. 

The plan of our paper is as follows. In section 2 we apply the optimal control 
theory to optical polarization in the weak field response. "A priori" pulse shape 
will not be conjectured. An exact and unique solution to a maximum coher- 
ence peaked at target time is derived for isolated I2 molecule in optically thin 
media. A time reverse (phase cojugate) resonance Raman scattering providing 
for a feedback to the vibrational spreading is represented as a basic principle 
of the quantum control. Controlling the pure and mixed states is considered. 
In section 3 we reconcile the quantum machinery with a reasoning appealing 
to the classical kinematics and discuss the phase modulated field of the scat- 
tered radiation associated with the classical action of the recoiling wave packet. 
Section 4 illustrates the general formulation by the numerical simulation. In 
section 5 we conclude and outline the future prospects. 

2. Mapping the Matter to a Light Field 

A rigorous mathematical treatment of molecular "magnetron" necessitates 
to solve the quantum equations dealing on equal footing with vibrational mo- 
tions and nonadiabatic electronic transitions. 

The adiabatic dynamics of a diatomic molecule is governed by Hamiltonian 
operators Hi, = Tkin + V& and H x = Tkin + V x . Here T^m is the operator of 
radial kinetic energy, where the only active coordinate is a distance between 
atoms, Vb and V x model the adiabatic potentials of the B and X states respec- 
tively. Having the fastest vibrational motion in diatomic molecule, we neglect 
the more slow molecular rotations (and the fine structure as well), since those 
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transients are separated by time scales . For heavy molecule as I2, the vibra- 
tional transient occurs on a subpicosecond time, while the rotational transient 
occurs on a longer time scale of 10 ps. Then, we can limit our consideration by 
the former, because the angular correlations happen past the vibrational ones. 
In fact, the rotational correction will only require a modification of the Hamil- 
tonians and averaging over initial rovibrational states. As justified in the works 
' 9 , the inclusion of the molecular rotations does not abandon the dynamic 
quantum control of vibrational wave packet towards a desired goal. 

Throughout this paper the Condon approximation for dipole transition mo- 
ment jj, is used, with fi being independent of internuclear separation in molecule. 
This assumption is valid for the weak field quantum control, when the molec- 
ular ground state is being well localized. The electronic states are coupled by 
the electric dipole operator V bx (t) = jl£(t), where the laser field is 

S(t) = E{t)e iul ^ h + E{t)e' iu)at/h 

and its slow envelop E(t) of allocated duration T p is to be found under a con- 
strain limiting the pulse energy, 

J = JcItE(t)E(t). (1) 


This laser field should prepare the molecular transition to a desired goal. Let 
optical dipole target be a linear superposition of weighted Dirac (5-like en- 
velopes shifted on t a , 

G{t)=Y J G a 8{t-t a ). (2) 

a 

The wisdom of that representation to the resonance coherence will be obvious 
further. 

The Schredinger equation for the wave functions ty b and ^ x reads as 

ihif b = (H b -u )V b + V bx {t)y x 

ihi> x = H x V x + vl(t)V b . (3) 

The empty B state and populated ground X state (at zero temperature) will be 
of use as our initial condition ^(O) = 0, ^(O) = ^ x ,o- Also, the rotating 
wave approximation is utilized in eq (3) to avoid fast oscillations of optical 

frequency of the B-X transition. Then, the resonance interaction is given by 

* 

the slow amplitude V bx (t) = fi E (t). 
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Our objective consists in controlling the polarization, 

oo 

D xb (t) = n(Vt\V b ) =fij drh x (r, t)* b (r, t), 

o 

which in its turn manifests in the molecular optical susceptibility. The first 
correction to the transition dipole moment D xb (t) in the weak field regime is 

T p 

D xb (t) = -iJ drS 00 (t -t)E (t), (t > T p ). (4) 

o 

To characterize the molecular transition, it is useful to extract a temporal struc- 
ture factor independing of the excitation field envelope E(t). For this goal we 
introduce a wave-packet correlation function as 

5oo(r) = ^(< |e-^/ ?l |^,o) (5) 

This formula describes the optical polarization caused by the <5-like pulse of 
electric field E s (t) = 5{t). Modulation of the B-X transitions is formed 
while the $ Ii(J replica propagates on the B curve. The bound wave packet os- 
cillates spreading between turning points and permits the molecule to radiate 
the Raman scattering signal. In frequency domain, the erratic spectral pattern 
17 of the resonance Raman intensity is accordingly observed. The spectral pro- 
file of the S correlator is represented by the KHD formula for the fundamental 
Raman overtone. Its resonance dependence is given by 



J 1 1 rr. 



7 + KPm - 



(6) 



rn 



where 7 is the damping constant; the eigenstates ^ bjm and eigenleveles j3. 
can be found from eigenequation H b ^ b m = l3 m ^ birn . The Stokes overtones 
are manifested when the wave packet reaches a favorable position to overlap 
the vibrational X states. 21 

To maximize the overlap between the Franck-Condon density l-D^i)! 2 and 
optical dipole target G(t), we define the field functional as 

Tp Tp 

F = J dtG(t)\D xb {t)\ 2 = Jdrj dnVoo^n) E(t) E(n). (7) 



One should note, that given functional does not confine the optical polarization 
throughout the time. However, it does guarantee a maximal dipole moment at 
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our objective G(t). The Frank-Condon transitions develop freely beyond the 
target time constrained only by the field energy in eq 1. The dynamic quantum 
control begins with the X state ^ x ,o, which is tightly localized in the I-I dis- 
tance at the start. The B state must be optimally driven below the dissociation 
limit to avoid bound-free transitions resulting in losses. Evident wave nature 
will impede the control, in which the B wave packet having undamped oscil- 
lation between turning points, must be focused on the inner steep core with 
a maximum velocity at target time. Those requirements will be met under a 
global maximum of the field functional F/J = X. This condition is set by 

small variation of the probe field E (t) in the variational equation, 

S(F - A J) = 0, (8) 

where A is the Langrange multiplier enabling to enforce the energy constrain. 
Herefrom the basic problem of dynamic quantum control in the weak field 
regime reduces to the Fredholm eigenequation, 

J dnVooir, n)E(ri) = \E(t), (9) 
o 

where the kernel of the homogeneous equation is given by the complex hermi- 
tian matrix, 

Poo(t,ti) = J dtG(t)S 00 (t - T)Soo{t - ri). (10) 

The integral kernel V depends exclusively on a chosen dipole target and speci- 
fication of the electronic transition being a pump-independent in the weak field 
limit. In some sense, it copies the material properties of the molecule to a light 
field. The time-dependent matrix elements eq (5) forming the kernel in eq (9) 
have long been known in the "matter-radiation" interaction theory, benefited to 
understanding nonlinear optics and utilized in quantum control of wave packets 
with different material objectives: the minimum space variance of wave pack- 
ets (I.S.Averbukh, M.Shapiro 3 ), the delta-like space density target (V.Dubov, 
H.Rabitz 7 ), the minimum quantum uncertainty state (K.R.Wilson et al. 5 ' 6 > 9 ). 

To recall their argumentation, we shall take a quick look at another field 
functional, which describes a total B population in the weak field regime, 

oo Tp Tp 

N = j dr\*£\r,t)\ 2 = J J drdriMoo(r,ri) E( Tl ) E{t), (11) 



where the self-conjugate kernel M.qq(t, t\) in eq 5 is proportional to the S- 
correlator depending on the difference argument t — t\, because the molecular 
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potentials are independent of time 

M 00 (t, n) = «S o(n - T)/h =s 00 (t - n)/n (12) 

Again one can repeat above steps searching a global maximum for the func- 
tional N/ J = X, which gives a maximum population yield per unit field en- 

ergy. Considering a small field variation E (t) for the variational equation 
5(N — A J) = 0, we obtain an eigenequation 

J dn S 00 (t - ri)^(n) = HXE(t) (13) 
o 

We obtain the degenerated kernel having separated the time arguments r and 
ri by means of the unity decomposition ^a,m^tt,m = 1- Then, a solution 

m 

to the integral equation can be cast as 

Eij) = -J* S(<ol**,m) e«W»^ (14) 

Unkno wn coefficients 5 m are found from the system of linear equations taken 
at a maximal eigenvalue X max 

A / lm,n£n — ^rnax&m (15) 

The matrix M. m n is formed by the Franck-Condon factors and energy levels 
of the B state 

e i{Pm-p n )T p /n _ i 

M m , n = (n,J^0)(*I, |* b ,n) i((3 m -p n ) (16) 

The optimal field E{t) depends on the oscillator strengths of the molecu- 
lar transition in eq(16), which includes factors of the increasing frequency due 
to the energy ladder /?o < Pi— < (3 m < ... of vibrational levels. This field 
must steadily populate the molecular states. It is expedient to put on the text- 
book Condon's model of the molecular transition between the parabolic curve 
V x = Mil 2 r 2 /2 and flat continuum Vb = const, where the S-correlator is the 
following: 



«Sbo(t) = (1 + iat/2)- 1 ' 2 = A (t) e 



•i0o(t)/2 



The phase profile (f>o(t) = arctg{£lt/2) is provided by the spreading wave 
packet in continuum which is projected to the ground vibrational state. On 
using a slope molecular curve VJ, = — fr models the repulsive Franck-Condon 
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region and recoiling wave packet. The phase correction <f>(t) = <fio(t) + 
t(ft) 2 / (12 Mh) is explicitly derived in the model dealing with the momen- 
tum representation of wave packet dynamics. The amplitude correction turns 
out to be of Gaussian type and like the Debye - Waller factor is given by the 
ratio of the recoil energy R = (ft) 2 / (2M) to vibrational quantum HQ, i.e. 

A(t) = A (t)e- R ^ 2hQ \ 

* 

The flat and slope continuum of molecular states provide for the kernel Soo(t— 
T\) with negative chirp of frequency for a short time duration t\, when the ex- 
pansion in the power of t\ is legitimate. The conclusion holds also in a gen- 
eral picture of semiclassical approximation as demonstrated in next section. 
Thus, a short pulse having positive frequency chirp slows down the fast inte- 
grand oscillations in eq (13) enabling one to maximize its eigenvalue, (i.e. the 
yield of the population). The maximum population has to involve a gain of 
fluorescence. 22 " 24 

Other trends of wave packet correlations make it possible to control a se- 
lected molecular state or coherent transients. According to eq 9, a top eigen- 
value A gives a maximum optical dipole transition per unit of field energy 
having eigenvector E\(t) as the globally optimal field of the allocated dura- 
tion T p . To demonstrate this statement we address the 5-Dirac model as the 
simplest optical dipole target, where the coherent envelope G(t) = S(t — td) 
will be peaked at time td after turning off the excitation pulse. The condi- 
tion td > T p means, that the "spontaneous" polarization will stand out on the 
pulse-free background. The delta-like target model results in the degenerate 
Fredholm kernel. An unique solution to the integral equation is explicitly writ- 
ten as 

E x (t) = S m (t d -T)(J/\) 1/2 , (17) 

where td marks a moment when the optical polarization is regained with A 
being the normalized coherent yield, 

A = J (It Soo (t d ~ T)S 00 (t d - t). (18) 
o 

The indication of optimality according to Bellman principle is indeed real- 
ized: the control depends on the state of system at the current moment alone. 
This globally optimal fields can be understood in general terms, for their en- 
velopes E\(t) match those of the wave packet correlations representing the 
time-reversed resonance scattering. The delay time td fits an absolute max- 
imum of the optical transient to target time r = td- The dipole transition 
moment is given by the convolution integral 
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D xb (t) = -i ( J/A) 1/2 J drS 00 (t - t) Soo (U - r). (19) 
o 

The optimal field in eqs 17, 19 cancels out the fast oscillating behaviour of 
the integrand as being designed to have the phase conjugate wave-packet cor- 
relations excited by the 5-like pulse. The optical transient by virtue of eq 19 
exhibits periodical recurrences, having the maximum magnitude — i(JA) 1 / 2 
exactly at the target time td- 

The straightforward extension of the 5-dipole target to the realistic shape 
G(t) specifies an eigenvalue problem in eq 9. The integral operator V can be 
decomposed into sum of the degenerate kernels according to eq 10 in the weak 
field limit. The solution is 

E ( T ) = J2 S oo(t a ~ T)G a £ a /X max (20) 

a 

where the £ a being the eigenstate at a top eigenvalue of the system 

Paai^ai = X max £ a (21) 

The matrix coefficients are given by, 

V aai = G a J dTS 00 (t a - t) Soo (t ai - t) (22) 
o 

so that the optimal field satisfies the energy constrain by the definition. 

We have established an important fact, which merits to be reformulated as a 
reversibility theorem asserting sufficient conditions to reverse the spread wave 
packet in time: In order to drive molecular transition to the polarization tar- 
get G(t), the optimal field should be made of superimposed phase conjugate 
fields of the resonance Raman scattering of ultrashort pulses of amplitudes 
G a £ a /X max at the delay time t a . 

The resonance scattering radiation itself may seem to be "erratic" in time. 
The wave-packet collapse and spreading and destructive interference are re- 
sponsible for the apparent noise due to molecular "disorder" on account of 
potential anharmonicity, uncommensurate frequencies, curve crossing etc. But 
phase conjugate feedback with respect to molecular correlations (resonant scat- 
tering) allows the wave packet to be recovered at a fixed space-time point at 
will. The optical coherence can be restored not only for the Raman fundamen- 
tal overtone. Our treatment might be readly extended to controlling the high 
overtones of molecular resonance scattering. 
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We term the field enabling to squeeze the optical coherence as a quantum 
hologram to note its key role for a wave packet interference in the Franck- 
Condon region. The quantum holography is relying on the wave packet corre- 
lations, which can be found by detecting a fluorescence or population created 
by two phase locked ultrashort pulses. Accordingly, in the limit of weak field, 
the excited state B is obtained from eq 3, in the first order to the resonance 
interaction as, 



(t) = -ie ie - itA »/H xfi (0) - ie 2 e-^- T ^ h ^ x , (T). (23) 



This superposition is obvious analog of an incident and object waves in optical 
holography. To produce the quantum interference of wave packets, the first 
laser pulse at t = must be followed by the delayed pulse at t = T. If the 
laser pulses are phase-locked and their areas are small with respect to ir, the 
population N of the state B contains the contribution 14 of the one photon 
two-pulse interference, i.e. 



The alternative pathways in course of the Franck-Condon transitions depend 
on the delay time T between pulses. This dependence is being just required 
to controlling the wave packet motion. The phase locked laser pulses must re- 
trieve synphase and quadrature components of the S correlator. The LIF signal 
1 , ionization channel 25 , or other means for the matter- wave interferometry 
may be employed. Then, the use of programmable optics 26,27 and algorithms 
based on the reversibility theorem fed into computer codes can be made to tai- 
lor the optimal field. This quantum holography program must be realized to 
rebuild the well localized replica of the ground state or to squeeze optical co- 
herence as called on even for unknown molecular curves as well as collisional 
or intramolecular dephasing and relaxation. 

The reversibility theorem still holds for mixed states inevitably residing in 
statistical systems, for which irreversible processes hamper the quantum con- 
trol of the wave packets. By considering the nondiagonal density matrix of the 
resonance transition p a p = (^t^), the brackets (...) symbolize an ensemble 
averaging in quantum kinetic theory. Here, the populations p xx = (\&t\& x ), 
Pbb = i^l^b) and polarization pi> x = (&1& X ) adhere to the kinetic equations: 



oc 




(24) 



o 



ih-£ = H,p -i(Tp) 



(25) 
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f Pbb, 


Pbx \ 


*- 


V Pxb, 


Pxx J 





where the rectangular brackets [...] denote a quantum commutator Hp — pH of 
the operators which are represented in the matrix notation as: 

f 66 ' Y b * ) (26) 
V xbl H xx J 

A phenomelogical damping matrix T designates the overall losses in irre- 
versible processes for collisions, spontaneous radiation, etc. Experiments in 
solids, liquids and gas cell conditions indicate an existences of coherent vibra- 
tional transients of I 2 to tens of picoseconds 14 . Thus, the quantum control 
theory within the matrix density formalism, in which the relaxation will not be 
scrutinized further, may be applied. 

Just as for pure states, we take care of the optical polarization p bx at target 
time. One can immediately maximize the overlap between the optical transi- 
tion pb x and its dipole target in eq 2. By letting the polarization operator p bx 
be a complex valued quantity, we define an overlap functional as 

F=-nfdt Im{G(t) Tr{ Pbx }} = J*' dr Re {E (t)V 00 (t)} (27) 

where another field functional Vqo(t)} is easily obtained from Eqs. (25), (26) 
to be 

Voo (t) = j~Jdt G(t)e~ r ^ Tr { e -*-^/» [p xx - p bb ] e ^-r)H x /n } 

(28) 

The functional F is indentified as a work that a given target G(t) would pro- 
duce on the Franck-Condon transition. This statement can be readily under- 
stood, since according to eq 25 the field functional F can be rewritten as the 
populations difference due to the target field "action" 

F=- J dt ImTr [G(t)p bx {t)\ = h J dt Tr [p b (t) - 'p x (t)] = h(N b - N x ) 

To maximize F, we apply the standard variational procedure of eq 8. The dy- 
namic quantum control in the weak field limit reduces to a problem of a linear 
mathematical programming. As a bonus of the approach, no integral equation 
needs to be solved at all, because, even for any target shape G(t), the varia- 
tional equation 

S(F - XJ)/S E{t) = 
results in an unique solution: 



E( T ) = A-^ooO-) 



(29) 
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By using the empty B state and populated state X, p^ = p xx = po, as the 
initial conditions, we set that 



p 00 ( T ) = !L dt G{t)e' v ^Tr{e- i{t - T)flb/h p G e i{t - T ^} (30) 



and the maximum coherent yield A becomes explicitly known from the energy 
constrain represented as 



By choosing the pure state po = ^ x ^ x ,o f° r the initial condition, we recover 
the main result of eq 5, as it must. Matrix density formulation has the ad- 
vantage of being applied to irreversible physical systems of gas and condensed 
matter. The predissociation 1 and caging phenomena can be considered as well. 
What is more, the observation and control of wave packets is in the progress 8 ' 9 
and includes rearrangement of subtle chemical structures in solids and liquids. 
Having postponed the feasible generalization, we shall discuss in next section 
the reversibility theorem in the semiclassical uniform approach. 

3. Semiclassical approximation 

The semiclassical limit of quantum dynamics is characterized by classical 
equations, that require to be supplemented by the quantum initial conditions to 
take full account of the singularity h = 0. The former is clear, only classically 
possible paths joining Rq with R\ contribute to the 5-correlator. These paths 
are invoked by Hamilton's principle SW = 0, where the classical action is 



If the action is known, it enables us to find a momentum Po = — dW(Ro, Ri,t)/ dR$ 
and density of states 8Po/dR\. Appearance of the derivative of classical tra- 
jectory in determining the density distribution in phase space is not occasional 
28 . The quantum transition amplitude connecting the points Rq and R\ de- 
pends on Plank 's constant h and can not be deduced from the classical postu- 
late. Instead, the quantum grounds dictate 



K(R ,Ri,t) = {Rq^-^'^Rx) = (-2i7rh5R 1 /5P y 1/2 exp(iW(R ,R 1 ,t)/h). 

(33) 



'j.Jortner has drawn our attention to curve crossing effects in experiments reported by V.Apkarian [8] et 
al. and M.Chergui [17] et al. Our theory can be extended for B-(B",a,a') couplings affecting the quantum 
control. The results is the entanglement with continuum molecular states shortening optical transients. 





(31) 




(32) 
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The probability density, \K,\ 2 = (2ttK)~ 1 5Pq/5Ri, stems from the total num- 
ber of atoms reached the point 5R\ started out from the point R and having the 
classical action W = W(Ro, Ri,t). The eq(33) does not give the quantum am- 
plitudes correctly for certain space-time points. Suffering from known draw- 
backs, it fails in the classically forbidden region, caustics and turning points. 
The story is old as the quantum theory itself. However, for our goals, the semi- 
classical approach does at least serve to a clue idea about the quantum control 

29 

The classical action W is inherited by the phase of electric field promoting 
resonant B-X transition. By this, the optimal laser field is related with the S- 
correlator which projects the wave packet B to the ground state X. The latter 
is centered at Rq = 2.7 A in the I-I distances with the variance AR = 0.05 A 
being much less than the Franck-Condon region extension itself. Hence, the 
mere closed paths passing the point Rq contribute to the optimal field. We can 
present its phase by the classical action eq 33 expanded for ultrashort pulse as 

W(Rq, Ro, t d -T)= W(R , Ro,t d ) + U b (R )r - M(R ) 2 t 3 /3. (34) 

in the vicinity of inner turning point Ro. The potential Ub(Ro) plays the role 
of frequency detuning of the B-X transition. Owing to zero velocity Ro = 0, 
its linear chirp will be negligible as compared with the nonlinear chirping. In- 
stead, there is a significant quadratic chirp — M(Ro) 2 r 2 t as adominant feature 
of acceleration R = —M~ 1 (dUb(R)/dR) due to the steep inner core. By 
travelling a distance Sr in vicinity of Rq, the wave packet detunes down the 
resonant frequency on Sr(dUb(R)/dR). The greater Sr = Rqt 2 the smaller 
the frequency of B-X transition. This gives the negative chirp, which qualita- 
tively is in agreement with the known reasoning 6 ' 29 . In practice, for a more 
longer pulse duration, this simple classical expression may overestimate the 
chirp rate, because the particle also travels far from the steep potential in a less 
repulsive Franck-Condon region. Reducing of frequency modulation, as well 
as linear chirping owing to a wave packet velocity is quite possible. 

A care must be taken to the turning points, where the semiclassical transition 
amplitude postulated in eq 33 diverges dR\ / dPo = 0. This is the case for our 
molecular magnetron design, in which the wave packet must be focused just 
in vicinity of the inner turning point. This severity can be circumvented by the 
uniform approach based on the Wigner representation of quantum operators. 
To find the S -correlator density, we must overlap the initial and propagated 
distributions in the phase space s = (r,p), 

4 f 

\S(-t)\ 2 = ^J drd PPo (r,p) po(R,P) (35) 
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where R = R(—t;s,0) and P = P(—t;s,0) being the current coordinate 
and momentum of a particle experiencing the force —dU\ ) {R)/dR. The clas- 
sical integrals of motion are fixed by the initial conditions R(0; s, 0) = r and 
P(0; s, 0) = p in the phase space. According to the Newton laws, the particle 
travels as, 

R = P/M P = -dU b (R)/dR, (36) 
and the Wigner density function obeys the transport equation: 

dpp(s,t) p dp (s,t) dU b (r)dp (s,t) _ 

dt M Or Or dp ' 1 ' 

Hence, the solution is po(R(0; s, t),P(0; s, t)) propagating along the trajectory 
passing the point s in time t with the distribution for the ground X state, which 
is defined as 

oo 

Po(r, P ) = J A e *P/»*t jo(r + |)^ o(r _ |) = ±re- F <P'\ (38) 

— oo 

where 

F(p, r) = ((r - R )/AR) 2 + (pAR/h) 2 . (39) 

The Wigner function of the ground state X is sharply peaked at the point 
qo = (Rq, 0) having Gaussian variances (AR, %/ AR) in the coordinates and 
momenta, respectively. The classical orbits passing this stretched ellipse make 
a major contribution to the 5-correlator amplitude. 2 Then, one can apply 
the fastest descent method to estimate the amplitude in eq 35. We expand 
the exponent functional F(R, P) to second order in the declines from the path 
So(t) = (R(0; qo, t),P(0; qo,t)) passing the equilibrium point qo of the ground 
X state. The family of orbits in its neighborhood contributes significantly. 
Thus, this topology determines basically a quadratic form 

F(R, P)=F° + F s ° (s - q ) + ^F s ° Sl ( S - q )( 8l - q ) + .., (40) 

The function F° = F(so(t)) is taken on the path so(t). A small stirring of 
its initial condition yields the four first and six second derivatives of R, P 
with respect to two-dimensional indices s = (r,p) or si = (r,p). Then, the 
quadratic form is settled by the twelve related functions: 

R(0; q , t), P(0; q o ,t),R s (0; q ,t),P s {0; q ,t),R s>Sl (0; go, t),P,, ai (0; q , t), 



2 It is worth noting, that a central role of the classical trajectories for calculating matrix elements of high 
excited states has been well known since the thirties, due to Landau seminal contribution to the chemical 
collisions theory 30 . 
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Of course, these functions adhere to Newton's equations and their derivatives. 
The useful vector and matrix notations 

(41) 

simplify the resulting Gaussian integral: 

\S(t)\ 2 oc h-^detF)- 1 * exp (-F° + ^VP^V^j . (42) 

The derivatives of the classical path with respect to its starting point in the 
space s represent the states density already mentioned in this section. Now we 
see, that a given method leads to the determinant in denomenator to set a den- 
sity of paths contributing to the quantum amplitude. In fact, this determinant is 
dictated by the uncertainty principle for the conjugate variables like coordinate 
and momentum. The merit of the average density valid for all paths, having 
the quantum spreading fixed by the initial distribution, is evident at the caustic 
and turning points. Here, the zeroth derivative Rp(t;s,0) = provides for 
the concentration of the classical paths, for which the Van Vleck's determinats 
diverge 28 ' 29 ' 31 . However, our determinantal relation still remains to be a good 
solution for the quantum amplitude. 

It is necessary to look more closely at the semiclassical recipe to gain insight 
about a frequency modulation of wave-packet correlations. The 5-correlator is 
written as a trace of the polarization density matrix, see also eq 30, 

5bo(t) = p 2 h~ 1 Tr{p(t)} oc J drdp p(s,t). 

s 

With initial condition p(0) = po, a solution formally satisfies to 

P (t) = e -itH b /h~ Qe itH x /h (43) 

The transport equation for the nondiagonal matrix element pb x (t) is immedi- 
ately obtained with accuracy fi 2 to be sure containing the description of quan- 
tum interference on the Franck-Condon transition: 

dp(s,t) p dp(s,t) dU d {r) dp(s,t) gg(r) (AA , 
dt M dr dr dp %h ^ S ' tj " U (44j 

The difference potential Ud{r) = Ub(r) — U x (r) defines corresponding reso- 
nant frequencies. An average potential is the half sum of the molecular curves 
U a (r) = (Ub{r) + U x {r))/2. The equations of motion are 

R = P/M P = -dU a {R)/dR (45) 
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So that, the 5-correlations are represented as an integral of local density with 
fast oscillating exponential prefactor in the space s given by the expression, 

t 

, -if drU d (R(r;s,t))/h 

Tr{p(t)} = drdp e ° Po(s(t)) (46) 

s 

At the start, the function po(s(0)) = po(R(0; s, 0), P(0; s, 0)) has a quite sim- 
ple shape of Gaussian distribution. Its contour lines look like the stretched el- 
lipses. Then, the map exhibits a complicated topology due to winding motion 
around a minimum of the average potential and its anharmonicity. Together 
with the fast oscillating functional for the difference potential of the B-X states, 
it presents a challenge to calculating the trace integral representation for 'topo- 
logical singularities' of the turning point structure which is common for the 
semiclassical approximations of wave mechanics 31 . 

The solution to our problem is written down in Langrange coordinates S = 
(R, P) = (R(0; s, t), P(0; s, t)). By virtue of Liouville theorem, the phase 
volume in the jacobian transformation from Euler to Lagrange variables is in- 
variant, whilst molecular potentials are independent of velocity. In the Lan- 
grange picture 

i 

, -if drU d (R(r;S,0))/h 

Tr{p(t)} = J dRdPe ° p {S), (47) 

s 

the basic integrand becomes a function of R, P 

p(t) = eM-i*(S,t)) p (S). 

To expand the path functional $ in the exponent to the second order in depar- 
tures from the point qo, we use a series 



t) = <D° + $0 (s - q ) + h° SySl (s - q )(sl - q ) + .., (48) 

where 

t 

$° = *(s ,t) = J dTU d (R(T;q ,0))/h 



and the path s (*) = (-R(0; qo, t), -P(0; qo, t)) obeys eq 45. By re-using the 
vector and matrix notation again, one can represent the final expression as a 
two-dimensional Gaussian integral, 
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S(t) = (fi/h) 2 (det ft) - ! exp - - Em -1 */^ , (49) 

where the vector U and matrix II are given by 

(Ai?)- 2 + f$ r ° r /2, \ 
*$ r ° p , (AR/h) 2 + i $p P /2 ) ■ 

(50) 

The phase modulation is set by the function Additional corrections in 

the exponent are responsible also for the decay amplitude giving generalized 
Debye-Waller factors for wave packet recoiling. The phase and amplitude for- 
mulae may be simplified due to peculiar features of molecular curves, short 
times evolution etc. This approach would be of interest to address Landau- 
Ziner transitions, and the spectral lines structure from impact center to static 
wings may be uniformly exploited. An analytical formulation of the globally 
optimal fields has some advantage of being used to model the molecular wave 
packet correlations. On propagating in a space the laser pulse might be tai- 
lormade to mimic the known optimal fields. With this aim the high frequency 
filtering in optical guides, nonlinear frequency dispersion, reflection from holo- 
graphic film, gratings etc. giving the right temporal and spectral trends may be 
adapted to perform the quantum control. 

As an example, we will discuss a chirping in static wing of spectral molec- 
ular line. Accessed by the ultrashort pulse, the wave packet develops in the 
average potential U a (R) between points connected by the direct path. For a 
small r (a femtosecond scale), these points are close together on the trajectory 

R( T ; R,P) = R + tP - (r 2 /2M) (dU a /dR) 

The difference potential in eq 46 can be expanded for the small departures from 
the equilibrium point go, and we may rewrite the frequency in $ as a power 
series for a small r 

U d (R(r;q ,0) = U d (R ) - (r 2 /2M) (dU a /dR )(dU d /dR ) = 
= U d (R ) - (t 2 /AM) (dU b /dR ) 2 . 

Since the condition dU x {R)/dR = (zero force) should be met in R = Rq 
(at a minimum of X curve), the result is the quadratic chirp lowering the pulse 
frequency with a correction coefficient of 0.25 to its classical rate in eq 34. 
The correction of this sort is not surprising for the classical - semiclassical - 
quantum correspondence. Moreover, considering more closely the phase and 
amplitude modulation of wave packet correlations resulted from Utl^W in 
eq 49, one can find additional corrections due to recoil from the slope molecu- 
lar curve. The first non- vanishing terms of the power r and r 3 in the exponent 
are an exact match the quantum formulae in the Condon model (see Sec. II 



c/ = (<i> r ,<i> p ); n = 
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after eq 16). It is worth stressing that the wave packet is found in the co- 
herent superposition of two resonance states, i.e. B and X, and their mixing 
will reduce the chirp rate on the classically gained formula in eq 34, because 
the molecular force acts only through the curve B (in recoiling from the inner 
curve), and does not in the equilibrium point of curve X. 

This argumentation must be revised further for returning trajectories on a 
long time duration. In the case, our uniform treatment warrants a more careful 
study that will be reported elsewhere. For completeness sake, the basic ideas of 
quantum control should be examined in a quantum simulation, which is favored 
to the semiclassical treatment, if one cannot limit by one classical path. Thus, 
in the next section, we shall rely on the quantum computation leaving aside a 
detailed numerical analysis of the classical equations. 

4. Numerical Simulations 

The quantum model eq 3 for the motion along one active I-I coordinate and 
electronic transition belongs to the polynomial class of numerical complexity 
of N linear equations, where N is the number of grid points. There is more that 
one recipe with which to perform their computations 32 . The QR diagonaliza- 
tion requires the 0(N 3 ) operations. In specific cases, the QR algorithm, as a 
stable spectral decomposition, can be used for the wave packet propagation. 
The time implicit symmetrical schemes of Gaussian elimination 33 necessitate 
O(N) operations. The symmetrized split operator method does 0(Nlog(N)) 
steps based on the fast Fourier transform (FFT) 34 . The FFT propagates the 
wave packet similar to Feynman's path integral and alternates the coordinate 
and momentum representations involving the thorough dynamical picture in 
hand. These methods are in good agreement. The practical value of FFTs or 
Gaussian elimination is that they enable to get rid of the matrix diagonalization. 

The absorbing boundary condition has been used in the numerical simu- 
lations. The imaginary negative optical potential discriminates the outgoing 
wave packet at outer edge of the numerical grid and has no effect on the 
bound states. This feature is visible in Figure 2, where the molecular dynam- 
ics following impact resonance interaction on B-X transition is shown. The 
ground state replica transforms into a wave-packet which is scattering above 
the molecular dissociation limit and binding below it. The scattering states, 
shown in the foreground of Figure 2, are involved due to impact photodissoci- 
ation (i.e., the 5-like pulse action). This branch is not mixed with the bound- 
bound transitions responsible for coherent recurrences, while the wave packet 
recoils against the molecular walls and oscillates spreading for the B state an- 
harmonicity. 

The resonance scattering radiation in the wake of ultrashort pulse proceeds 
in two stages: a fast monotonic decay due to direct photodissociation gives 
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way to a "non-regular" reviving. This picture is shown in Figure 3 for the 
fundamental overtone of Raman intensity. Since a damping out is assumed to 
be small, the wave packet correlations hold during many molecular periods. 
The amplitude and phase profile will follow in reverse order. Thus, Figure 3 
showing its fundamental overtone intensity must be read " from right to left". 
The phase conjugate scattering signal in Figure 4(a) has been designed to drive 
the spread waves towards td- The target time td exactly represents the starting 
point of wave-packet correlations. Furthermore, on going back in time their 
envelope S(td — t) is not delta correlating as a white noise signal. 

There are two regular trends: the negative frequency chirp, which has been 
explained with the semiclassical arguments, and amplitude growing to the 
pulse end. The last feature stemmed from the reversibility theorem can be also 
understood qualitatively. Since the wave packet is spreading between turn- 
ing points, the result is a small amplitude almost everywhere in the classically 
accessible Franck-Condon region. The optimal field on its start is set by a pro- 
jection of the spread wave packet to the X state. The S -correlations undergo 
fast oscillations due to multiple recoils of the wave packet against molecular 
walls and its spreading accompanied by decreasing amplitude. Evolving back 
in time, the wave packet becomes more regular, and its shape is recovered and 
copies well localized replica of the ground X state. Focussing the wave packet 
in the vicinity of the inner turning point results in a further rise of its amplitude 
3 . The projection to the ground X state grows together. This explains why the 
globally optimal field behaves in such a way in Figure 4(a). Induced resonance 
transient given by eq 4 is shown in Figure 4(b). The 1 ps square gate of optimal 
field creates the transient spike at delay time 1.5 ps. The temporal full width at 
half maximum (FWHM) equals T w = 25 fs that is in good agreement with the 
kinematical overlap lifetime r rec = 15/s. 

This semiclassical estimate is astonishingly accurate although less then quan- 
tum computed. There are two main reasons of this distinction. The wave 
packet interference is absent in the classical kinematics and the dispersion of 
excited state is not taken into account. Thus, the ignorance of quantum wave 
nature underestimates the width of the coherent peak. To cure the situation, the 
semiclassical uniform approach, which copes with initial quantum distribution, 
classical motions, caustic and turning points, may be applied. 

For definiteness sake, setting T g = 50 fs square gate of optical target G at 
the delay time 1.5 ps, we now can consider a realistic optical transient. Av- 
eraging the globally optimal envelope over G washes out the high frequency 
components as shown in Figure 5. However, the average pulse shape is yet 



3 A clear analogy is given by a point of intense radiation behind burning-glass exposed to sun. 
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endowed with a chirp. This feature is exhibited in Figure 6 (a), where Wigner 
spectrogram is shown. The Wigner function of the optimal fields is defined as 



Picosecond square gate envelope and periodic boundary condition have been 
used in evaluating the spectrograms. The frequency-time maps manifest a 
striking correlation between the optimal fields and wave packet evolution. For 
instance, the envelop of allocated duration T p = Ips consists of more shorten 
subpulses. Their number is associated with the number of vibrational cycles 
between turning points during T p . For every recoiling from inner molecu- 
lar core the wave packet gives rise to a subpulse. The contour levels of the 
two-dimensional Wigner spectrogram are taken at one-half a maximum height. 
A slope and concavity characterizes respectively a linear and quadratic fre- 
quency chirp of subpulse. The semiclassical reasoning relate the chirping to 
wave packet velocity and acceleration. The chirp is in fact linear because 
time duration of pulse T p (in picosecond range) is much more than the re- 
coil lifetime r rec = 15fs. Thus, the wave-packet is far from turning points 
and does not accelerate most of the time. The spectral FWHM in Figure 6 
(a) is about Au « 500 cm -1 being consistent with the reciprocity relation 
Au = 2irh/AT g 35 . 

A minor decoherence of the objective as shown in Figure 6b is obtained in 
our numerical simulation. The temporal FWHM is about 60 fs being 10 fs over 
the square gate target width T g = 50 fs. The distinction stems from the wave 
packets quantum dispersion. The wave packet dynamics is displayed in Figure 
7 and Figure 8. The resonance frequency is 19050 crrT 1 above the ground X 
term; the excitation bandwidth is about 500 cm -1 (FWHM). Thus, the wave 
packet B is created below the molecular continuum to spread from 2.2A to 4A 
not far from its dissociation limit. Given intensity of transitions, the depopula- 
tion of state X is involved to the pulse end, although the X state localization is 
not affected at the saturation. The optical transient rises and falls, following the 
wave packet recoils from the inner core. The spatial variance of a maximum 
focused wave packet at target time reaches 0.04 A, being slightly less than the 
original localization in agreement with our qualitative reasoning. The result is 
robust to slow sinusoidal or Gaussian aberrations of optimal fields and akin to 
known stability of phase conjugate fields due to the wave speckle structure 36 . 

5. Conclusion and Future Perspectives 

The numerous possible applications in femtosecond spectroscopy are rested 
on the concept of coherence. The proceedings of Femtochemistry III Confer- 
ence in Lund testify (see the special issue of the Journal of Physical Chem- 
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istry), that this realm of knowledge is flourishing and another lines of re- 
searches have now been opened. As an example we refer the field of deemed 
quantum computers 37-39 . Here, controllable coherence of logic gates could 
boost performance of massive parallel computation exponentially with num- 
ber of quantum nodes. These nodes may be constituted from clusters or single 
molecules and their resonance optical transitions must act in concert with vi- 
brational wave packets captured by "molecular cavities". A run stage of a quan- 
tum computer, its coherence may not be affected by environment 40 . Thus, the 
danger of dephasing requires the optical coherence to be actively controlled ac- 
cording to a state-mapping protocol. From this point of view, the reversibility 
theorem stated in our paper exemplifies a paradigm of controllability expected 
for wave packets dynamics and optical transients. In fact, the fundamental 
connection between physical laws and computations based on time reversable 
Hamiltonians has been addressed yet by R.Feynman 41 . The computational 
quantum networks requires the reversable transmission among distant nodes 

42 

The central aim of this paper is to set a globally optimal solution to a max- 
imum squeezing the optical coherence of a single molecule. The sufficient 
conditions of recovering coherence have been stated in the reversibility theo- 
rem. Our time reversed fields have been tested in the model computation of 
I2 B-X transition confirming effeciency of the phase conjugate resonant Ra- 
man scattering for quantum control. The result is not limited only to diatomic 
molecules in gas or condensed states, but is valid for the quantum control of 
electron wave packets as well 43 > 44 . 

The four-wave mixing interaction and related stimulated Raman emission 
is deserving of mention. These optical transients would make it possible to 
automate the quantum control to be held both for bound states and molecular 
continuum. To take one example, ultrashort transform limited pulses driv- 
ing bound-free transitions result in the photon echo of photodissociation 45 , to 
mean an automated controlling wave packet motion. Photon echoes 19 imply 
that dynamical holograms are written on resonance Franck-Condon transitions. 
But it takes more strong laser fields to induce at least the third order nonlin- 
ear polarization. A self-adaptive "quantum mirror" has been devised 36,46 to 
reverse wave evolution on curing the optical aberrations of a noise wave front 
propagated trough inhomogenious media. Moreover we have seen that the 
phase cojugate (time-reverse) Raman scattering negates also derealization of 
molecular wave packets. 

The use of the phase conjugate signals presents an intriguing illustration 
of how the original quantum state is restored in the matter-wave interferome- 
try. An uniform semiclassical formula has been derived for the optimal field 
envelop. The simple estimates are contrasted with quantum simulation. We 
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conclude that femtosecond coherence of a single molecule must exist in varied 
conditions (even on resonant transitions saturated to the pulse end). Robustness 
of our molecular magnetron design has been demonstrated also. We propose 
to test experimentally the molecular magnetron design verifying the quantum 
holography directly on observing delayed sparks of femtosecond fluorescence 
which is excited by a more long picosecond pulse in optically thin sample. 
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Figure 1. The resonance curves of h molecule [20]. Our design represents the iodine mag- 
netron. A tailored laser pulse, as shown by the point-filled box, excites the vibrational wave 
packet. Then, the wave packet is reflected from outer curve and moves back to focus on the 
inner Franck-Condon region. The result is a tight overlap with the ground state X. The maximal 
optical polarization must be followed by the resonance fluorescence burst at target time. 
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Figure 2. The picture shows the nodal structure of bound wave packet ranging between 
2.2 — AAA. The wave packet fringes are due to interference of the wave packet in virtue of 
multiple collisions against the molecular curves. 
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Figure 3. Intensity of fundamental overtone of resonant scattering as a function of time. This 
resonance transition is induced by the photon impact (the delta-like pulse). The scattering B 
states contribute to monotonic free decay signal first few fs. The bound states interference, 
present in the background in Figure 2, is responsible for reviving the resonance Raman scatter- 
ing after 320 fs delay, which is matched with a molecular period (at 570 nm in B state). 
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Figure 4. The square 1 ps gate of the optimal field in (a). This optimal field creates the dipole 
transition moment shown in (b). The accomplishment of our coherent objective is evident at 

target time, pmifll to 1 S 
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Figure 5. The intensity of optimal field in (a) and its real part in (b). This optical pulse is 
designed to create 50 fs square gate optical coherence at target time 1.5 ps. 
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Figure 6. The frequency-time plot of the globally optimal pulse for the 50 fs square gate 
target is given by the contour map at half a maximum of Wigner spectrogram in (a). The 
coherent transient in (b) is created by the optimal field. The temporal FWHM equal to 60 fs fits 
well the 50 fs square gate target. 
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Figure 8. The wave packet B free propagates after the pulse to a tightly squeezed state at the 
inner turning point to overlap well the ground state X at target time 1.5 ps. 



